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ABSTRACT 

We analyze the physical properties of post-starburst AGN 
SDSS J085338. 27+033246.1 according to its optical spectrum and discuss 
its implications on AGN's evolution. The spectra PCA method is developed to 
extract emission lines and absorption features from the total light spectrum. 
The emission-line analysis indicates that the object can be classified as a NLS1 
with FeII/H/3 B = 2.4 ± 0.2, large Eddington ratio (~ 0.34), small black hole 
mass (~ 1.1 x 10 7 M Q ) and intermediately strong radio emission. A simple 
SSP model indicates that the absorption features are reproduced by a ~100 
Myr old starburst with a mass of ~ 7 x 1O 9 M rather well. The current SFR 
~ 3.OM yr _1 inferred from the [Oil] emission is much smaller than the past 
average SFR~ 7OM yr" 1 , however. The line ratio diagnosis using the BPT 
diagrams indicates that the narrow emission lines are almost entirely emitted 
from HII regions. We further discuss a possible evolutionary path that links 
AGN and starburst phenomena. 

Subject headings: galaxies: active — galaxies: starburst — quasars: emission lines - 
quasars: individual (SDSS J085338. 27+033246.1) 
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1. Introduction 

In recent years, accumulating evidence supports the hypothesis that AGN reflects an 
important stage of galaxy formation. Observations reveals the fact that AGN activity and 
star formation frequently occur together (e.g., Cid Fernandes et al. 2001; Kauffmann et 
al. 2003; see reviews in Gonzalez Delgado 2002). An important clue of co-evolution of 
AGN and starburst is the fundamental, tight relationship between mass of supermassive 
black hole (SMBH) and the velocity dispersion of bulge where the SMBH resides in (e.g., 
Magorrian et al. 1998; Gebhardt et al. 2000; Tremaine et al. 2002, Greene & Ho 2006). 
By analyzing the spectra of narrow-line AGNs from Sloan Digital Sky Survey (SDSS, York 
et al. 2000), Heckman et al. (2004) found that most of accretion-driven growth of SMBH 
takes place in the galaxy with relatively young stellar population. Recently, as essential 
relationships relating to AGN's activity, the Eigenvector 1 space (Boroson & Green 1992, 
hereafter BG92) was found to be related with ages of the stellar populations as assessed by 
infrared color a (60,25) (Wang et al. 2006). 

It is then natural to ask a question "What is a young AGN?" . Mathur (2000) argued 
that Narrow-line Seyfert 1 galaxy (NLS1) with high Eddington ratio and small black hole 
mass might be a young AGN. If so, the co-evolution of AGN and starburst implies that 
NLSls are expected to be associated with relatively young stellar populations. 

To understand the elusive AGN-starburst connection, we need to carefully examine 
the properties of AGN's host galaxy. However, optical spectroscopic study is possible for 
only a few type I AGNs because the strong nuclear emission usually masks the relatively 
faint stellar components. By fitting deep off-nuclear optical spectra Canalizo & Stockton 
(2001) and Nolan et al. (2001) determined the ages of circumnuclear stellar populations 
for a few IR and optical-selected QSOs respectively. In the happy case, Brotherton et al. 
(1999) discovered a spectacular "post-starburst quasar" UNJ1025-0040 with a ~400 Myr 
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old stellar population. Wang et al. (2004) identified SDSS J022119.84+005628.4 as a NLS1 
(RFe=1.4±1.0) whose spectrum displays a clear Balmer jump. However, the authors did 
not perform further stellar population analysis. 

In this paper, we perform a detailed examination of the physical properties of 
post-starburst AGN SDSS J085338. 27+033246.1 and discuss its implications on study 
of AGN's evolution. The optical spectrum is extracted from SDSS Data Release 4 
(Adelman-McCarthy et al. 2006). Its appropriate spectrum allows us to separate it into 
a young stellar population and a typical NLS1 spectrum. In addition, the emission line 
profiles allows us to determine the location on the line ratio diagnostic diagrams. We 
begin with a brief description of the data reduction in §2. The results and implications 
are presented in §3. The cosmological parameters: h = 0.7, Qm = 0.3, and Q\ = 0.7 are 
adopted for calculations. 



2. Spectral Analysis 

We extract the optical spectrum of SDSS J085338. 27+033246.1 when we carry out a 
systematic search for broad line AGNs from the objects listed in the catalogues 1 released by 
Kauffmann et al. (2003) and Heckman et al. (2004). Totally hundreds of such objects are 
extracted from the parent sample. The details of the search and properties of the sub-sample 
will be presented in subsequent papers. Among these objects SDSS J085338. 27+033246.1 
has unique properties with strong enough typical AGN spectrum superposed on a spectrum 
of A-type star. It is coincident that we find the object has been identified as a member of a 
sample of 74 post-starburst broad line AGNs by Zhou et al. (2005) when our paper is being 
prepared. 

1 These catalogues can be downloaded from http://www.mpa-garching.mpg.de/SDSS/. 
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The observed spectrum is smoothed with a boxcar of 3 pixels (~5A) to enhance S/N 
ratio and precision of the spectra measurements. Standard IRAF procedures are adopted to 
reduce the raw data, including Galactic extinction correction (with E[B — V) = 0.048 from 
NASA/IPAC Extragalactic Database (NED) and assuming R v = 3.1) and de-redshift to 
rest frame with z = 0.207881. The total light spectrum at rest frame is displayed in Figure 
1 (the second line from top to bottom). The spectrum clearly shows a prominent Balmer 
jump at blue even diluted by an AGN continuum and high order Balmer series. 

In the next step, the principal component analysis (PCA) method is developed to 
remove the stellar light from the observed spectrum (e.g., Li et al. 2005; Hao et al 2005). 
We build a library of stellar absorption spectra by applying PCA technique on standard 
SSP models developed by Bruzual & Chariot (2003, hereafter BC03). The first seven 
eigenspectra are used to model the starlight component by their linear combination. In 
order to appropriate model the starlight component, our template contains the seven 
eigenspectra, a power-law continuum, a Fell complex of AGN (adopted from BG92) and a 
Galactic extinction curve (Cardelli et al. 1989). In order to avoid the distortion of low S/N, 
a x 2 minimizing is performed over the rest wavelength range from A3700 to A6800, except 
for the regions around the strong emission lines. The removal of the stellar component is 
illustrated in Figure 1 as well. Also plotted are the modeled Fell blends and emission-line 
spectrum of the AGN. The flux of the Fell blends measured between rest frame wavelength 
A4434 and A4684 is (4.23 ± 0.10) x 10~ 15 ergs s" 1 cm" 2 . The detailed comparison between 
the observed and modeled spectrum at blue is shown in the insert plot A. It is clear that 
the modeled spectrum reproduces the observed absorption features rather well. 

The AGN emission lines are modeled by SPECFIT task as described in Wang et al. 
(2005, 2006). Briefly, each narrow line (e.g., Ha, H/3, [Nil]) is modeled by a Gaussian 
component. The intensity ratio of [Nil] doublets is forced to equal to the theoretical 
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prediction. The spectra modeling is schemed in the insert plot B and C for H.(3 and Ha region, 
respectively. The upper limit of [OIII]A5007 flux FQOIII]) < 1.2 x 1CT 16 ergs s" 1 cm~ 2 
is estimated by integration rather than fit because the feature is very faint. In 
addition to these lines, the spectrum displays a broad emission feature at about 5880A. 
Based on its asymmetric profile, we identified this feature as HeIA5876 contaminated 
by broad NaIDAA5890, 5896 emission. We find EW(HeI+NaI D)~ 5.6 ± 1.3A and 
(HeI+NaID)/Ha B = 0.11 ± 0.02. The NalD emission was detected in a few AGNs (e.g., 
Veron-Cetty et al. 2004; Veron-Cetty et al. 2006 and references therein, Thompson 1991). 
The calculations indicated that the existence of large NalD emission (NaID/Ho;=0.01-0.05) 
is related to models with high density (N e ~ 10 11 cm~ 3 ) and large column density 
(A H > 10 23 ' 5 cm~ 2 ) (Thompson 1991). The measured line properties are summarized in 
Table 1. The propagation of error is included in the uncertainties shown in parentheses. All 
the uncertainties given in Columns (2) and (3) are caused by profile modeling. 



3. Results and Discussion 

3.1. Emission Line Analysis: A NLS1 with Intermediate Strong Radio 

Emission 

According to its spectra properties, we classify SDSS J085338. 27+033246.1 as a NLS1 
both because of the extremely faint [OIII] emission ([OIII]/H/3 <0.05) and because of the 
relatively narrow B./3 broad component (FWHM(H/3b) = 1827.6 ± 113.8 km s _1 ). Moreover, 
the measured RFe is as large as 2.4 ± 0.2, which is typical of NLSls. 

We further calculate its black hole virial mass and Eddington ratio in terms of the Ha 
broad component according to the equations in Greene & Ho (2005): 

BH V10 42 er S s s- 1 / V 1000 km s- 1 / ° KJ 
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L 5100 = 2.4 x 10 43 ( in42 LHa r ) °' 86 erg (2) 
\10 42 erg s 1 / 

The intrinsic luminosity L(Ha B ) = (2.59 ± 1.02) x 10 42 ergs s" 1 is used for calculations 
because of the large intrinsic extinction E(B — V) = 0.42 ± 0.18. The extinction is 
calculated from the line ratio Hon/H/Sn = 4.8 ± 1.0, assuming a Galactic extinction 
curve and Ry = 3.1 (Osterbrock 1989). The inferred virial mass and Eddington ratio is 
~ 1.1 x 1O 7 M and ~ 0.34, respectively. It is clear that all these properties (i.e., small 
black hole mass, large RFe value and large Eddington ratio) indicate the object should be 
potentially classified as a NLS1. 

The radio radiation at 1.4GHz of the object was detected by NVSS (Condon et al. 1998) 
and FIRST Survey (Becker et al. 1995). The map of the FIRST survey shows an unresolved 
source with a flux ~4.04mJy and position a = 08 h 53 m 38 s .269, 5 = +03 ° 32'46".39 (J2000). 
The discrepancy in position given by the optical and radio surveys is much less than 1" 
(~0".3). We determined the extinction-corrected continuum flux that is contributed by 
AGN at rest wavelength A4400 to be 0.16 ± 0.02mJy. If the radio emission is mainly due to 
AGN, the inferred radio loudness R defined as R = F iadio / F X uoo (Kellermann et al. 1989) is 
logi?i. 4G Hz=l-40 corresponding to logi?=1.12 for 5GHz radio flux 2 . Adopting the division 
\ogR = 1 between radio loud and quiet AGN, SDSS J085338.27+033246.1 is an object with 
intermediately strong radio emission. 

We further calculate its radio luminosity Pi.4gh z — (5.0 ± 0.3) x 10 23 W Hz -1 . 
This luminosity has same order of magnitude of the most radio-luminous starburst 
(logP 4 .85GH z — 22.3 — 23.4, Smith et al. 1998). It has long been known that the decimeter 
radio luminosity emitted from star-forming galaxy traces supernova rate of massive enough 
stars (i.e., M > 8M ). The lifetime of these massive supernova progenitors is ~ 10 7 ' 5 yr. The 

2 -Ri.4GHz = l-9i?5GHz assuming a spectral shape f u oc z/~ 5 from optical to radio band. 
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radio supernova itself and its remnant has lifetime ~ lOOyr and ~2x 10 4 yr, respectively. 
The average SFR over the past ~ 10 7 - 5 ~ 8 yr therefore relates with the radio emission 
attributed to supernova as (Condon 1992) 

SFR <^> = 4.0x^WH Z - M ° yr " < 3 » 

For the object, the estimated average SFR is ~ 70 M & yr _1 (see Sect. 3.2 for details). 
Using a Salpeter IMF (Salpeter 1955), this corresponds to a SFR(> 5M Q ) ~ 14M yr -1 . 
The inferred luminosity is therefore ~ 5.6 x 10 22 W Hz -1 and is about one order lower than 
the total radio luminosity. 



3.2. Post-starburst: "A-type star" spectrum and residual star formation 

The "A-type star" spectrum indicates the age of the stellar population is no more than 
a fewxlO 8 yrs. To interpret the spectrum of the young stellar population, we extract a 
series of spectra of SSP models of BC03 with a solar metallicity and with a Chabrier initial 
mass function. The extracted spectra are matched with the observed spectrum 3 one by one 
within the wavelength range 3500-5000A. The spectrum red of A5000 is not used because 
this part is obviously contaminated by an old population. Figure 2 illustrates the best 
match between the observed spectrum and a O.lGyr old SSP model. The model predicted 
mass of the starburst is ~ 7.0 x 10 9 M Q . This simple SSP model shows a past average SFR 
~ 70 M yr -1 . 

However, recent studies suggested that the AGN and starburst activity perhaps not 
co-evolve simultaneously. Schmitt (2001) suggested that the starbursts may predominate 

3 In this study the spectrum used for age diagnosis is corrected for intrinsic extinction 
derived from the PCA modeling. 
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over AGNs in their earlier phase. Yan et al. (2005) recently suggested that most of SDSS 
post-starburst galaxies hold somewhat AGN like activities. The post-starburst is the phase 
in which the galaxy already ceased star formation at some recent epoch. The derived 
~ O.lGyr post-starburst in SDSS J085338. 27+033246.1 is somewhat similar with the 
~ 0.4Gyr post-starburst that was obtained by Brotherton et al. (1999) in UN J1025-0040. 

[OII]A3727 line is widely used as an empirical indicator of ongoing SFR (< 10 7 yr) for 
emission line galaxy surveys (e.g., Gallagher et al. 1989; Kennicutt 1992; Hippelein et al. 
2003). Its practicability led to several calibrations and a great deal discussions (e.g., Kewley 
et al. 2004; Kennicutt 1998; Yan et al. 2005). The extinction-corrected flux of [OII]A3727 
doublets is (1.9 ± 0.7) x 10~ 15 ergs s _1 cm -2 (with its observed flux shown in Table 1), 
which yields a ratio log([OII]/[OIII])>0.55. Adopting the ratio log([OIII]/H/3 N )<-0.65, 
the object is consequently located in the region that represents star formation (Kim et al. 
2006 and references therein). It means that the [Oil] emission ought to be almost entirely 
attributed to star formation (see below). 

We use the calibration 

SFR - 7 ' 9X 16.73 -1.7 L So/H) + 12] < 4 > 

(Kewley et al. 2004) to estimate the SFR, where £[oii],42 is the luminosity of [Oil] emission 
in units of 10 42 ergs s -1 . The SFR is inferred to be round about 3.0 M yr _1 by adopting 
the assumption used in Ho (2005), i.e., the metallicity is twice of solar corresponding to 
log(0/H) + 12 = 9.2. Here, we obtain a residual but quenched ongoing SFR as obtained in 
the optically selected QSOs by Ho (2005). The distance- independent birthrate parameter b 
is defined as the ratio of the current SFR to the average past one (Kennicutt et al. 1994). 
As a particular case, the object has extremely small value of b ~ 0.04, i.e., has present SFR 
significantly lower than the past. 



3.3. Location on the BPT diagram 



The traditional BPT diagram (Baldwin et al. 1981) is a powerful tool for diagnosing 
the origin of emission of narrow lines for emission-line galaxies. The location of the object 
in the BPT diagram in which the line ratio [OIII]/H/3 is plotted against [Nil] /Ha is marked 
by solid square in Figure 3 (left panel). The solid and dashed line shows the demarcation 
line between AGN and starburst galaxy defined by Kewley et al. (2004) and Kauffmann 
et al. (2003), respectively. The figure shows clearly that the position of the object is 
marginally above the threshold used to define representative HII region. The similar plot is 
shown in Figure 3 (right panel) but for [OIII]/H/3 v.s. ([SII]A6719+6731)/Ha. In this case, 
the position of the object falls into the region that represents starburst and is far below 
the demarcation line. The combination of these two diagnosis suggests the observed narrow 
emission lines are mainly produced by star-formation. What is the distribution in the BPT 
diagram for post-starburst AGNs is a further interesting question. 



3.4. Implications 

Here it is interesting to give more discussions on the elusive AGN-starburst connection. 
Wang et al. (2006) recently indicated that the well-documented El sequence most likely 
represents an evolutionary track of AGN. The track further implies that AGN with high 
Eddington ratio evolves to one with low Eddington ratio. This result also naturally implies 
that NLSls that occupy one extreme end of El should be associated with relatively young 
stellar populations. Apart from the evidence mentioned above, for the IR-selected QSOs, 
Canalizo & Stockton (2001) found the ages of circumnuclear stellar populations range from 
current star formation to ~ 300Myr old. Subsequent spectroscopy studies indicated that 
about 70%-100% IR-selected QSOs are extremely or intermediately strong Fell emitters 
(e.g., Lipari et al. 2003; Zheng et al. 2002). Moreover, Zhou et al. (2005) found that more 
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than half of the post-starburst type I AGN fulfill the criterion for NLS1. All these evidence 
appear to reveal that a young AGN is most likely a strong Fell emitter with relatively 
narrow H/3 profile. In contrast, old stellar populations (~ 8 — 14Gyr) are found to dominate 
the off-nuclear (~ 5") stellar population for optical-selected QSOs (Nolan et al. 2001 and 
references there in). 

These accumulating clues lead us to propose a possible evolutionary scenario that links 
both AGN and starburst phenomena. The massive starburst might be dominant in the 
earlier evolutionary phase. Then the starburst passively evolves to A-type stars when the 
plenty gas falls into the center of a galaxy under the gravitational attraction of SMBH. 
At one stage about ~ a few x 10 8 yr after the beginning of the starburst, the SMBH begin 
to accrete matter with small black hole mass and high Eddington ratio. The feedback 
of the AGN suppresses circumnuclear star formation at the same time (e.g., Kim et al. 
2006). The black hole mass grow substantially in the formed bulge in this short phase 
(e.g. Mathur et al. 2001; Mathur & Grupe 2005). The duration of the phase is commonly 
estimated as e-folding time scale t = 4.5 x 10 7 (^-)(-j-^-)yr (Salpeter 1964), where i] is the 
radiative efficiency The differential growth of AGN and starburst allows the NLSls below 
the Mbh — o"* relation for inactive galaxies to approach normal Mbh — o"* relation (Mathur 
& Grupe 2005; Mathur et al. 2001; Wandel 2002; Bian & Zhao 2004). 

The AGN then shines with decreasing accretion rate and insignificantly increasing 
black hole mass because of the consumption of the gas. The scenario implies that the 
AGN and host co-evolve marginally around the common Mbh — < J * relation. The luminous 
AGN phase persists for typical AGN lifetime which is usually believed to be 10 7 ~ 8 yr with 
however large uncertainties (see reviews in Matini 2004). During this phase, the young 
stellar population continuously ages and its importance gradually fades. The underlying 
old stellar population (a few Gyr) finally becomes dominant in the emission of host of an 
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old AGN. In fact, in addition to the old stellar populations, relatively old post-starbursts 
(0.1-2Gyr) are frequently detected in off-nuclear regions of a few powerful radio- loud AGNs 
(e.g. Tadhunter et al. 2005 and references therein). After the luminous phase, AGN then 
passively evolves to a stage with low luminosity, or low radiative efficiency 

In summary, this scenario could be described as following path: starburst — > 
NLSl+post-starburst — > luminous AGN with decreased Eddington ratio+old stellar 
population — > less luminous AGN. 

We believe future studies of AGN+post-starburst composite objects such as 
SDSS J085338. 27+033246.1 are important for developing a rational scenario to interpret the 
AGN-starburst connection. 



4. Conclusions 

The physical properties of post-starburst AGN SDSS J085338. 27+033246.1 are derived 
by analyzing its optical spectrum. It allows us to obtain results as following. The object 
can be identified as a NLS1 (([OIII]/H/3 < 0.05, FWHM(H/3) = 1827.6 ± 113.8 km s" 1 and 
RFe=2.4±0.2) associated with a post-starburst stellar population as identified from the size 
of the Balmer jump. A simple SSP model indicates the starburst with a mass 7 x 1O 9 M 
possibly took place ~ 1.0 x 10 8 yr ago. The [Oil] inferred SFR is as small as ~3.OM yr _1 
which is much smaller than the past average one ~ 7OM yr~ 1 . Its locations in the BPT 
diagrams show the HII region significantly contributes to the emission of the narrow lines. 
We further discuss its implications on the elusive AGN-starburst connection. 
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Fig. 1. — Subtracting the stellar component using the linear combination of the 7 eigenspec- 
tra and modeling of the emission line profiles. From top to bottom we plot the emission line 
spectrum, observed spectrum overlaid by the modeled spectrum, modeled spectrum of the 
host galaxy and modeled Fell blends. The insert plot A shows detailed comparison between 
the observed and model spectrum at blue. The emission line fit is displayed in the insert 
plot B and C for H/? and Ha region, respectively. 
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Fig. 2. — An illustration of the match between the stellar spectrum and a O.lGyr old SSP 
model with solar metallicity and Chabrier initial mass function. The SSP model is shifted 
downward arbitrary for visibility. 
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Fig. 3.— Left panel: The location of SDSS J085338.27+033246.1 on the BPT diagram 
defined by line ratio [OIII]/H/3 v.s. [Nil] /Ha. The density contours are shown for typical 
distribution of the narrow-line galaxies described in Heckman et al. (2004) and Kauffmann 
et al. (2003). Only the galaxies with S/N > 20 and the emission lines detected at at least 3<r 
significance are plotted. The solid line is the demarcation line used to separate AGN from 
starburst by Kewley et al. (2004), and the dashed line by Kauffmann et al. (2003). Right 
panel: The same BPT diagram but for [OIII]/H/3 plotted against ([SII]A6716+6731)/Ha. 



- 21 - 



Table 1: Properties of the emission lines of SDSS J085338.27+033246.1. Flux of each 

component is normalized to the flux of H/3 narrow component F(H/3 N ) = (5.4 ± 0.8) x 

10~ 16 ergs s _1 cm -2 . The uncertainties including propagation of errors are shown in the 
parentheses. 



Line identification 


Flux ratio 


FWHM 






(km s _1 ) 



H/? N 1.0 301.8 ±38.0 

H/fe 3.2±0.2(±0.6) ... 1827.6 ±113.8 

[OIII]A5007 <0.22 

Ha N 4.8±0.4(±1.4) ... 314.9 ±23.3 

Ha B 16.1 ± 1.1(±3.1).. 1778.6 ±150.0 

[NII]A6583 2.7±0.3(±0.6) ... 402.2 ± 36.8 

Fell 4570 7.8 ± 0.2(±1.4) 

[OII]A3727 0.56±0.02(±0.10) 

[SII]A6716 0.59±0.06(±0.12) 240.2 ± 10.0 

[SII]A6731 0.35±0.01(±0.06) 

HehfNalD 1.85 ± 1.24(±0.41) 



